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Abstract 

Thermal measurements on a GaAs/AlGaAs heterostructure reveal that 
the state of the confined two-dimensional electrons dramatically affects the 
nuclear-spin diffusion near Landau level filling factor i^=l. The experiments 
provide quantitative evidence that the sharp peak in the temperature depen- 
dence of heat capacity near z^=l is due to an enhanced nuclear-spin diffusion 
from the GaAs quantum wells into the AlGaAs barriers. We discuss the phys- 
ical origin of this enhancement in terms the possible Skyrme solid-liquid phase 
transition. 

PACS numbers: 73.20.Dx, 73.40.Hm, 65.40.-hg 
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In the presence of a perpendicular magnetic field (B) , two-dimensional electron systems 
(2DESs) exhibit striking phenomena originating from the Landau quantization and electron- 
electron interaction; examples are the integral and fractional quantum Hall effects At 
Landau level filling factor z/=l, theoretical work on 2DESs in GaAs/AlGaAs heterostruc- 
tures has shown that electron spin textures known as Skyrmions are the lowest energy, 
charged excitations of the ferromagnetic ground state 0]. Skyrmions, which result from the 
dominance of the electron-electron exchange energy over the Zeeman energy, are also the 
lowest energy state for quasi- holes and quasi-electrons near v=l\ the ground state of a 2DES 
near v=l is in fact expected to be a crystal of Skyrmions 

Recent experimental studies have revealed that Skyrmions play a major role in many 
properties of quantum Hall systems 0-0. In particular, measurements on a GaAs/AlGaAs 
multiple-quantum-well (QW) heterostructure indicated an anomalous behavior for the low- 
temperature apparent-heat-capacity (C) near z/=l: C is enhanced by a factor of up to ~10^ 
with respect to its low-i? value 0. For temperatures T > O.IK, the magnitude as well as 
T- and 5-dependence of the data are consistent with the large C being dominated by the 
Schottky nuclear heat capacity of Ga and As atoms in the QWs. This observation implies a 
strong enhancement of the nuclear spin-lattice relaxation rate near z/=l, consistent with the 
results of nuclear magnetic resonance (NMR) experiments and recent calculations |^ 
which attribute the enhancement to the presence of low-energy spin-flip excitations in the 
2DES. An even more striking feature of the C vs T data is a very sharp peak at a very low 
temperature Tc(z/), suggestive of a phase transition in the 2DES 0. While recent estimates 
for the melting temperature of the Skyrme crystal 0,0, including its dependence on u, 
suggest that the peak in C could originate from a Skyrme solid-liquid phase transition, the 
physical mechanism that affects the nuclear spin dynamics and gives rise to the anomalous 
C-peak is not known and is a matter of debate 

In this Letter we report new quasi-adiabatic thermal experiments revealing that the 
mechanism responsible for the peak in C vs T is a dramatic enhancement of nuclear spin 
diffusion across the QW-barrier interface. While only the nuclear heat capacity of the QWs' 
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atoms is observed away from Tc, the nuclei in both the QWs and barriers contribute to the 
measured heat capacity near T^. We discuss the physical origin of this phenomenon in terms 
of the possible Skyrme solid-liquid phase transition. 

The experiments were performed on the same multiple-QW heterostructure used in 
Ref. p. It consists of 100 GaAs QWs separated by Alo.3Gao.7As barriers. The wells and 
barriers are 250 and 1850A thick, respectively, and the barriers are 6 doped with donors (Si) 
near their centers. We used a 7 x 7mm^ piece of the wafer which was thinned to 65/im. The 
lattice T is probed by a carbon paint resistor deposited on the substrate side of the sample 
and electrically connected with four very low thermal conductance, 8/im-diameter, NbTi 
wires. This resistor was calibrated against a RUO2 resistance thermometer at B=0. Since 
the nuclear spin temperature (T/v) may be inhomogeneous in the sample due to a variable 
nuclear spin-lattice relaxation rate , T strictly refers to the measured lattice temperature. 
The sample is supported by two 60/im-diameter NbTi wire which serves as the main weak 
thermal link to the cold finger of the dilution refrigerator mixing chamber. In the restricted 
T and u-range investigated here, C is very large and the observed external time constant. 
Text, reaches extremely large values which exceed 10"^ seconds. Therefore, in the time scale 
of our experiments, which is governed by the observed internal time constant Tint ~ 10^ s 
(see below), the sample is in the quasi-adiabatic regime. B is applied at an angle of 6' = 30° 
with respect to the 2DES plane in order to match the experimental conditions of Ref. 0. 

We first summarize the heat capacity results of Ref. . Figure |l] shows the measured 
C vs T at z/=0.77, revealing a sharp peak at Tc=42mK. While at high T (^70mK) the 
measured value and T-dependence of C are rather consistent with the calculated Schottky 
nuclear heat capacity of Ga and As atoms in the QWs P,p!0| (solid curve in Fig. [l|), at lower 
T, and especially near T^, C exceeds the calculated value by a factor of up to ~10. The 
peak value of C, however, appears consistent with the Schottky nuclear heat capacity of the 
heterostructure if the nuclei of the ftarner atoms are also included (dotted curve in Fig. 
This observation suggests that, near Tc, the barriers' nuclear spins also contribute to the 
measured C. 
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In our new quasi-adiabatic experiments, we utilize another noteworthy feature of the 
heat capacity data; as seen in Fig. |I] inset, Tc sensitively depends on filling factor. This 
dependence, which is reproduced in Fig. ^ (a) (dashed curve), is the key to understanding 
our results. In these experiments, we first fix the magnetic field at B=8.5T (z/=0.77) where 
the peak in C occurs at Tc=42mK (Fig. |T] and dashed curve in Fig. 0(a)). Starting from 
T~60mK, we cool down the cold finger to base temperature (~10mK) and the lattice slowly 
cools down to T~32mK by waiting for about 50 hours. Then, in order to further improve the 
adiabatic conditions, the weak heat flow that results from the T-difference between the lattice 
and the cold finger is reduced by increasing the cold finger temperature to 32±3mK. Under 
these conditions, we observe that the measured sample T is stable within the experimental 
accuracy (^ ±0.2mK) over periods longer than Text, meaning that the sample is effectively 
quasi-adiabatic. 

In a first experiment, whose results are shown by open circles in Fig. @(a), we started from 
such a quasi-adiabatic condition and swept B from 8.5 to 7.7T, at a rate of 0.01 T/minute. 
We observe that T rises from 31.7mK to 36.2mK. B was then swept back to 8.5T at the 
same rate; this led to a rise of T to 39.5mK. While the increase in T with increasing B may 
be explained by the adiabatic magnetization of the nuclei in the QWs, which couple to the 



lattice, according to |10 



Tn oc B, (1) 

the T rise observed when B is lowered from 8.5T to 7.7T is puzzling. 

In order to better understand the results of the above experiment, we performed a similar 
experiment, starting from T=32.6mK, but here we swept B in steps of 0.2T from 8.5T to 
7.3T, with a hold time of 3 hours between steps. The evolution of T during this second 
experiment is presented in Fig. ^(b) and summarized in Fig. |^(a) by close circles. We note 
that in each step, except for steps 2^3 and 3—^4, the lattice T rises during the -B-sweeps 
and then decays to approximately the same value as the one at the end of the previous 
step. These T rises are due to Eddy-current heating of the 2DES and the lattice during 



the -B-sweeps. The heat is evidently slowly absorbed by the nuclei during and following the 
sweep via internal relaxation. Because of the very large heat capacity of the nuclear spins, 
this (Eddy-current) heating does not lead to an appreciable increase in the lattice T at the 
end of the step |Tl|]. The key feature of this experiment, however, is that during the 2— s-3 
and 3— i>4 steps, the lattice T increases significantly while the sample is in quasi-adiabatic 
conditions. Note in particular that, during the 2^3 step, the lattice T rises even after the 
i?-sweep is completed. This striking behavior occurs when the lattice T approaches Tc at 
which C exhibits a maximum (dashed curve in Fig. 1(a)). Indeed, it appears that near Tc 
the lattice is heated internally. 

To explain the results of Figs. || and 0, we propose a model in which the barriers' nuclei 
couple to the QWs' nuclei very near T^, but decouple at higher and lower temperatures. 
The sharpness of the peak in C vs T data (Fig. P is consistent with this model and implies 
that the coupling between the QWs and barriers is turned off rapidly as T deviates from 
Tc. Therefore, we expect that during the initial cool-down at B=8.5T, when T decreases 
below Tc=42mK, T/v in the AlGaAs barriers remains close to 42mK while T/v of the QWs 
decreases down to f»32mK due to strong coupling to the lattice. Next, decreasing B below 
8.5T (steps 1^7) in adiabatic conditions has two distinct consequences: (1) T/v in both the 
QWs and barriers is reduced due to adiabatic demagnetization of the nuclei (Eq. P, and (2) 
the coupling between the QWs and the barriers will increase dramatically when the dashed 
(Tc vs B) curve in Fig. H(a) is crossed. While the former is essentially monotonic in B, the 
latter is not and implies that, near the dashed curve in Fig. 0(a), T/v should equalize over the 
whole GaAs/AlGaAs heterostructure. Since the barriers' nuclei were initially warmer than 
the QWs' nuclei, we expect a rise in lattice T near this crossing, as observed experimentally. 
Finally, we attribute the T rise during the final (7— >8) step to the nuclear magnetization of 
the nuclei (Eq. |I|). 

Beyond the good qualitative description of the experiments in Fig. 0, the above model 
appears to provide a reasonable quantitative account of the data also. First, when B is 
lowered from 8.5T to 7.3T, we expect from Eq. |l] that nuclear demagnetization reduces T/v 
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in the barriers from :^42mK to ?5i36mK. Since the barriers are by far thicker than the QWs, 
according to our model, this value should approximately correspond to the temperature of 
the entire heterostructure, and hence of the lattice. This prediction is in good agreement 
with experiments that give T^37mK at B=7.3T. The final increase of T from Ri37mK to 
?s43mK when B increases from 7.3T to 8.5T is consistent with the nuclear magnetization of 
the heterostructure (Eq. Second, the fact that, at Tc, C reaches the value expected for 
the entire heterostructure (dotted curve in Fig. |l|) is clearly consistent with our model and 
hence provides additional credence to our interpretation. 

Next, we discuss the internal time constant (Tint) observed in our experiments. Figure |^ 
presents a typical heat pulse/relaxation trace, in quasi- adiabatic conditions, obtained during 
heat capacity experiments near T^. The relaxation follows an exponential decay characterized 
by Tint- The inset to Fig. ^ shows that Tint exhibits two different, but remarkably constant, 
values above and below 30mK: ~1600s and ~900s, respectively. We note that this crossover 
temperature corresponds to the T below which C is reduced back to a value close to the 
nuclear heat capacity of the 100-QWs, Cqw (full curve in Fig. |I]). 

In order to explain these observations, we consider three possible mechanisms responsible 
for the observed Tint- Heat diffusion in the nuclear spin system of the barriers is governed 
by the one- dimensional diffusion equation: 



where D ^ 10 ^''m^/s is the nuclear-spin diffusion constant in GaAs [T^ and the z-axis is 



along the growth direction. Therefore, it gives an internal time constant: 



„2 



~ ^ (3) 



where r is the distance over which diffusion takes place |]T3|. Since each barrier is surrounded 
by two QWs, we take r=925A (half the barrier thickness) and Eq. ^ gives r^i^SOOs. The 
two other mechanisms involved in the heat transfer from the lattice to the barrier nuclei, i.e. 
nuclear spin-lattice relaxation in the QWs and diffusion across the QW-barrier interfaces, re- 
sult in two additional time constants: Ti and Tj, respectively. When C ~ Cqw {T ^ 30mK), 
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the diffusion across the QW-barrier interface is neghgible and Tmt should be determined by 
Ti only. On the other hand, when C > Cqw (30 ^ T ^ 50mK) at least a fraction of the 
barriers' nuclei couple to the QWs' nuclei. This implies that and should increase Tint- 
This is consistent with the rise in Tint observed above 30mK (Fig. ^ inset). Moreover, we 
note that the increase of Tint above 30mK (~700s) is comparable to the estimated rd~800s, 
implying that Ti is rather small once the barrier nuclei do couple to those in the QWs 



We now remark on the physical origin of the peak in C vs T . NMR experiments show 
that nuclear spin-diffusion from the QWs into the barriers is extremely weak, except when 
optical pumping broadens the Knight shift peak of the QWs which then overlaps with the 
NMR peak of the barriers 0. The spectral overlap allows spin-diffusion which is driven by 
nuclear magnetic dipole-dipole coupling [|6|,p!5|. Therefore, the enhancement of spin-diffusion 
across the QW-barrier interface near could originate from a broadening of the Knight 
shift peak in the QWs. In the liquid Skyrme phase, motional averaging of the local 2DES 
spin polarization produces a single Knight shift peak [^0- On the other hand, the absence 
of motional averaging in the solid Skyrme phase should induce both positive and negative 
Knight shifts depending on the local spin polarization of the 2DES. Since above and below 
the peak in C vs T the diffusion across the QW-barrier interface is very weak, this implies 
that in both the liquid and the solid Skyrme states, the Knight shift peak(s) do not overlap 
with the NMR peak of the barrier. One possible explanation for the C anomaly is that 
the critical slowing down of the spin fluctuations in the 2DES, associated with the Skyrme 
liquid-solid phase transition 0, could induce a broadening of the Knight shift peak in the 
QWs. This broadening would induce spectral overlap of the QWs' and barriers' NMR peaks 
and hence allow spin-diffusion across the QW-barrier interface only near T^. 

Finally, it is worth emphasizing that the shape of the peak in C vs T points to a second- 
order phase transition |T6|. The fact that Tint is constant in the T-range where C > Cqw 
near Tc implies that the strength of the coupling to the barriers remains constant while C 
varies strongly. Therefore, it appears that the smaller value of C on the sides of the peak in 
C vs T comes from a reduced fraction of the sample in which the barriers' nuclei couple to 



the QWs' nuclei. The complexity of spin textures in the 2DES, together with the fact that 
the phase transition is not observed directly but through its effect on nuclear-spin lattice 
relaxation, does not allow a quantitative analysis of the data in the framework of the existing 
theories that describe order-disorder phase transitions. On the other hand, Skyrmions form 
an XY ferromagnet which has broken translation symmetry and one expects at least one 
phase transition as T is lowered 

In conclusion, our experiments reveal the critical influence of the ground and excited 
states of the 2DES near z/=l on nuclear-spin dynamics in GaAs/AlGaAs heterostructures. 
Our results provide a quantitative phenomenological description of the apparent heat ca- 
pacity anomaly that may originate from a Skyrme solid-liquid phase transition. 
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950328 and the NSF MRSEC grant DMR-9400362. V.B. acknowledges financial support of 
the Belgian National Fund for Scientific Research. 
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FIGURES 

FIG. 1. Measured C vs T at i?=8.5T and 9 = 30° {u = 0.77). The curves represent the 
calculated Schottky nuclear heat capacity of the GaAs 100-QWs (full curve) and the 100-period 
GaAs/Alo.sGao.rAs heterostructure (dotted curve). The inset shows the measured z^-dependence 
of Tc at 6* = 30°. 

FIG. 2. Measured lattice T during the quasi-adiabatic experiments, (a) summarizes the evolu- 
tion of T during the first (o) and the second (•) experiments (see text). The full and the dotted 
lines are guides to the eye and the arrows indicate the direction of evolution. The dashed curve 
corresponds to the S-dependence of as reproduced from Fig. |l| inset, (b) shows T and B vs 
time in the second experiment. -B-sweeps (gray areas) separate hold times at the end of which T 
is recorded (1 to 7) and summarized in (a) together with the final T at i?=8.5T (8). 

FIG. 3. T vs time during a heat capacity experiment in the quasi-adiabatic regime [C = Q/AT 
where Q is the heat injected during the pulse; i?=8.5T and 9 = 30°). The heat pulse is followed 
by a relaxation characterized by an exponential decay (full curve) giving Tint ~ 1500s. The inset 
shows Tint vs T at the same B and 9. 
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